On the extreme positive star-formation feedback condition in 

SCUBA sources 



Sergiy Silich^, Guillermo Tenorio-Tagle^'^, Casiana Munoz-Tunon'^, Filiberto 
Hueyotl-Zahuantitla^, Richard Wiinsch^ 

and 

Jan Palous^ 

ABSTRACT 



We present a detailed study of the hydrodynamics of the matter reinserted by massive 
stars via stellar winds and supernovae explosions in young assembling galaxies. We show that 
the interplay between the thermalization of the kinetic energy provided by massive stars, 
radiative cooling of the thermalized plasma and the gravitational pull of the host galaxy, lead 
to three different hydrodynamic regimes. These are: a) The quasi-adiabatic supergalactic 
winds, b) The bimodal flows, with mass accumulation in the central zones and gas expulsion 
from the outer zones of the assembling galaxy, c) The gravitationally bound regime, for 
which all of the gas returned by massive stars remains bound to the host galaxy and is likely 
to be reprocessed into futher generations of stars. Which of the three possible solutions 
takes place, depends on the mass of the star forming region its mechanical luminosity (or 
star formation rate) and its size. The model predicts that massive assembling galaxies with 
large star formation rates similar to those detected in SCUBA sources (~ 1000 Mq yr~^) 
are likely to evolve in a positive star-formation feedback condition, either in the bimodal, 
or in the gravitationally bound regime. This implies that star formation in these sources 
may have little impact on the intergalactic medium and result instead into a fast interstellar 
matter enrichment, as observed in high redshift quasars. 
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1. Introduction 

In Cosmology today the study of the star forma tion negative feedback is recognized as 



one of the central issues regard i ng galaxy formatioii ( iDekel &: Silklll986t iFriaca fc Terlevich 



19981 : IScannapieco et al.l 120021 : iFerreras et al.l |2002[ ) . The large UV photon output from 



massive stars and their violently deposited mechanical energy, make them indeed major 
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Tenorio-Taele & Bodenheimer 


shown by 


Tenorio-Tagle et al. 



Wiinsch et al.l (120081 ) for the case of massive and compact super stellar 
clusters, this may not be the whole story as the stellar feedback may, in extreme cases, 
become positive. This would allow gravity to win over thermal pressure. This may also 
be the case if one considers star formation in the submillimeter (SCUBA) galaxies - high 
redshift sources with the highest (~ 1000 yr~^) star formation rates (SFRs) so far 
known. These have been observed in the submillimetre continuum (emission from warm 
dust in the rest-frame far-IR/submillimetre wavelengths), in the CO line emission associated 
with the cold molecular gas and in the near-infrared integral field spectroscopy, which deals 
with the rest-frame optical emission lines asso c iated with the pho t oionized gas ( see, for 
example. 



Hughes et al. 


1998; 


Greve et al. 


2005: 



and references therein). 

The global properties of SCUBA galaxies include a typical dynamical masses of 5 ± 
3 X 10^^ M0, a gas to dynamical mass fraction, fg = Mgas/Mdyn ~ 0.25 — 0.3, and a 
radius of about 1-3 kpc, par ameters consistent with the expected properties of massive 



spheroids in the early Universe (iGreve et al.ll2005t ISwinbank et al.ll2006t iTacconi et al.ll2006 



Schinnerer et al.ll2008[ ). The long star formation duty-cycle with a time-scale ~ 100—300 Myr 
and the inhomogeneou s nature of SCUBA s ources favor a continuous star formation scenario 
( ISwinbank et al.l l2006t iTacconi et al.ll2008l ). Here we show how massive and violent star- 
forming events driven by a high rate of star-formation lead to positive feedback. Looking 
at some extreme cases, we identify radiative cooling as the agent capable of downgrading 
the impact of the stellar energy deposition, leading inevitably to an extreme positive star 
formation feedback condition which should play a major role in galaxy formation. We also 
show that gravity is another major player. The gravitational pull of the galaxy also leads to a 
positive feedback condition, particularly in compact proto-galactic sources. The gravitational 
pull then prevents the formation of a supergalactic wind, retaining the injected and ablated 
matter within the star forming region, favoring its accumulation and conversion into further 
generations of stars. 



In section 2 we examine the physical implications of massive star formation rates and 
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discuss the physical hmits between the various possible hydrodynamic regimes. The impli- 
cations of events with a high SFR and our conclusions are given in sections 3 and 4. 



2. Star formation under a large SFR 



If one scales the evolutionary synthesis models (e.g. iLeitherer fc Heckmanlll995l ) for star 
clusters generated by a const ant star forniation rate to the values inferred from the SCUBA 
sources (> 100 Mq yr'^; e.g. ICreve et al.lboosi ISwinbank et~aDl2006[ IXacconi et aUbooeh . 
one sees that as a result of the continuous death and creation of massive stars, the UV 
photon output will level off at ~ 10^^ ionizing photons s~^ after 3 Myr of the evolution. 
The mechanical energy deposited by the evolving stars (Lsf) through winds and supernovae 
(SNe) will also increase, although not so rapidly, to reach a constant value ~ 2.5 x 10^^ erg 
s~^ after 40 Myr of evolution. Accordingly, the mass violently returned to the ISM by stellar 
winds and supernovae will amount to 3 x 10^ Mq after 10 Myr, reaching almost 10^ Mq after 
100 Myr of evolution. The absolute values of all the above mentioned variables ought to be 
linearly scaled by more than an order of magnitude, at the given times, if instead of a SFR 
equal to 100 M0 yr~^, one assumes the even larger values inferred for the most powerful 
SCUBA sources (> 1000 Mq yr'^). 

At first glance, such an energy deposition and such a vast amount of matter so violently 
injected, would u navoidably lead to extreme massive outflows into the intergalactic mediurn 
see, for example, Heckman et al. 199ol: Strickland fc Stevens 2000; Scannapieco et al. 2002 



Tenorio-Tagle et al.l 



2003 



Veilleux et al. 



2005 



and references therein). Supergalactic winds 
are believed to result from the full thermalization of the kinetic en ergy of the ejecta, throug h 
multiple random collisions within the star-forming volume (see IChevalier fc Cleggi Il985l ). 
Thermalization generates the large overpressure that continuously accelerates the deposited 
matter to flnally blow it out of the star-forming volume, composing a stationary super- 
wind with an adiabatic terminal speed Vaoo = C^Lsf/ MsfY^"^', where Lsf and Msf are 
the mechanical energy and mass deposition rates provided by stellar winds and supernovae 
explosions within the star-forming volume. For this to happen, the ejecta has to reach 
an outward velocity equal to the sound speed (c = {'jP/pY^'^ oc T^/^) right at the star- 
forming boundary, Rsf, to then fulflU the stationary condition in which the rate at which 
matter is deposited equals the rate at which it streams away from the star-forming region: 
Msf = 4:hR'^pPsfCsf, where psF and csf are the values of density and sound speed at the 
surface of the star forming region. However, as shown in the present series of papers (see 
Sihch et al.l[2003l . booi iTenorio-Tade et al.lEoOsl . boO?! : Iwiinsch et allboosh . when dealing 
with the outflows generated by massive bursts of star formation, the impact of radiative 
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cooling becomes a relevant property, as is gravity, able to hold a fraction of the deposited 
matter within the star cluster volume. 

In the case of an instantaneous burst of star formation, stellar winds and supernovae 
are able to remove the matter left over from star formation out of the star cluster volume 



in just a few Myr (IMelioli fc de Gouveia dal Pinol l2006l : iTenorio-Tagle et al.l |2006| ) and so 
the hydrodynamic solution considers only the matter reinserted by the massive stars. In 
the continuous star formation scenario however, a gas reservoir out of which a constant SFR 
could be sustained is required. This implies that besides the mass returned by supernovae and 
stellar winds, Msf, the flow may hold additional matter. This results from the destruction 
and mass ablation from star forming clouds and can be normalized to the star formation rate 
within the star forming region: M/^ = rjidSFR , where riid is the mass loading coefficient. 
The total mass input rate into the flow is then: 



M = Msf + Mi 



Id 



2Lo 

^Aoo 



+ riid X IMeyr" 



SFR 
IMoyr- 



'1^ 



where Lq is the normalization coefficient, which relates the mechanical energy output rate, 
LsF, with the SFR: 

LsF = L^{SFR/lMQyi"^). (2) 

Hereafter we shall adopt Lq = 2.5 x 10"^^ erg s~^ and Vaoo = 2750 km s~^ These values result 
from Starburst 99 synthetic models for a continuous star formation mode with a Salpete r 
IMF and sources between 0.1 and 100 M©, for ages t > 40 Myr (ILeitherer et al.lll999l ). 
Note, that mass loading changes the outflow terminal speed, which in this case is smaller 
than Vaoo- 

2LsF \ ^^'^ Vaoo 



Msf + Ml 



Id 



1 + 



2Lo 



1/2- 



(3) 



For the calculations our sem i-analytic radiative stat ionary wind cod e (ISilich et al.ll2008 ) 



and th e ID hydrocode ZEUS3D (IStone fc Normanlll992l ) as modifled by iTenorio-Tagle et al. 



(120071 ) and IWiinsch et al.l ( 120081 ) were used. All numerical calculations were performed with 
an open inner and outer boundary conditions. The hydrodynamic equations here include the 
gravitational pull from the baryonic and dark matter, both assumed to be homogeneously 
distributed inside the star forming region. In our approach it is also assumed that the mass of 
the flow is negligible compared to the dynamical mass of the system, and thus the self-gravity 
of the re-inserted gas is not included in the calculations. For example, the mass of the flow 
within the star forming region, Mjiow, normalized to the dynamical mass of the system is: 
Mfiow/Mdyn ~ 2 X 10"^ and Mfiow/Mdyn ~ 8 x 10"^ in case of model 1 with a low SFR and 
model 4 with a high SFR, respectively. To obtain the stationary hydrodynamic solution one 
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has to know the mechanical energy and mass deposition rates, which are defined by equations 
([1]) and ([2]), and the radius of th e star formin g region, Rsf- We use the equihbrium coohng 
function, A(T,Z), tabulated by iPlewal (Il995l ) and set the metallicity of the plasma to the 
solar value in all calculations. Our reference models are presented in Table 1. Here the first 
column marks the model in our list, the ablation coefficient, rjid, is presented in column 2, 
columns 3, 4 and 5 present the radius, dynamical mass of the star forming region and the 
star formation rate, respectively. Column 6 provides information regarding the resultant 
hydrodynamic regime. 



Table 1: Reference models 



Model 


md 


Radius 


Dynamical mass 


SFR 


Regime 






(kpc) 


(10" Mq) 


(Mq yr-i) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1 


0.5 


2.5 


2 


2 


Superwind 


2 


0.5 


1.65 


2 


2 


Superwind 


3 


0.5 


1.2 


2 


2 


Gravitationally bound 


4 


0.5 


2.5 


2 


1200 


Bimodal 


5 


0.5 


1.7 


2 


1200 


Bimodal 


6 


0.5 


1.2 


2 


1200 


Gravitationally bound 



3. The hydrodynamic regimes 

There are three major hydrodynamic regimes that develop within galaxies undergoing a 
large SFR. Which of the three possible solutions takes place, depends on the mass of the star 
forming region and its position in the mechanical luminosity or SFR versus size (Rsf) pa- 
rameter space. Figure 1 presents the threshold lines, which separate proto-galaxies evolving 
in different hydrodynamic regimes. The left panel presents threshold lines for proto-galaxies 
whose dynamical mass is equal to 2 x 10^^ Mq for different values of the ablation coefficient 
ivid = 0.1,05,1.0, dotted, solid and dashed lines, respectively). Below the threshold lines 
radiative cooling has a negligible effect on the flow and the reinserted matter ends up as a 
superwind. Above these lines radiative cooling leads to a bimodal regime in which some of 
the reinserted matter within the densest central regions loses its pressure and is unable to 
parti cipate in the galactic wind. Instead it accu r nulate s there fueling further stellar genera- 
tions |TenOTioJagk_et.alJ (120051 ): IWiinsch et al.l (120081 ) . For compact star forming regions. 



to the left of the vertical lines shown in Figure 1 [Rsf < Rcrit), gravity inhibits the formation 
of a superwind, leading instead to matter accumulation and to further generations of star 
formation. In these cases the sound speed at the surface of the star-forming region is smaller 
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than ~ [GMdyn/^RspY^'^ , which is one half of the escape speed from the proto-galaxy sur- 
face. The threshold lines for less (4 x 10^° Mq) and more (4 x lO^-*^ M©) massive galaxies with 
?7;d = 0.5 are presented in Figure 1, right hand panel. 




Fig. 1. — The threshold SFR or energy input rate versus the threshold size. The left panel 
displays the threshold mechanical luminosity, SFR and critical radii for different values of rjid 
in the case when the dynamical mass in the star forming system equals to M^yn — 2x 10^^ Mq. 
The limiting energy input rate and its corresponding constant star formation rate (right-hand 
axis) above which strong radiative cooling inhibits the stationary superwind solution, as a 
function of the size of the star-forming region and the mass ablation coefficient, rjia. These 
are terminated at the vertical lines, which display the critical radii, Rcrit- Gravity inhibits 
the formation of supergalactic winds in systems with smaller radii. The right hand panel 
shows how the location of the threshold lines depends on the total mass of the star forming 
region. Several of the cases here presented are marked by crosses in the left panel. Note that 
the vertical line that marks the critical radius (Rcrit) for rjid ~ 0.1 in the left panel, and for 
Mdyn — 4: X 10^° Mq, in the right panel, lies within 1 kpc and thus is not shown. 



As initially expected (see section 2) there is a large fraction of the parameter space that 
leads to stationary supersonic winds. In these cases, all the deposited matter, as well as that 
ablated from clouds, is able to escape from the gravitational well of the galaxy. For this to 
happen the flow has its stagnation point (the point where the velocity of the flow equals to 
km s^^) right at the center of the galaxy (Rgt = pc) and its sonic point at the surface. 
The matter accelerates then through pressure gradients to reach supersonic velocities and 
form a supergalactic wind as it streams away from the galaxy. 

Models 1 and 2 in Table 1 undergo such supergalactic winds. The distribution of the 
hydrodynamical variables in this case is shown in Figure 2. Here panel (a) presents the 
flow velocity (solid line) in the case of model 1 and compares this to the local sound speed. 
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('jP(r) / p{r)Y^'^ and Vesci^) = (2GM(r)/r)^/^, which is the escape velocity at any distance 
r outside the star-forming region. The outflow velocity reaches the local sound speed value 
right at the surface of the star forming region, then it accelerates rapidly to reach its terminal 
value of ~ 740 km s~^ at a distance about 4 kpc from the galaxy center. At this distance it 
already exceeds the escape velocity and thus composes a supergalactic wind. Panels (b) and 
(c) present the distributions of temperature and density in the flow. The temperature drops 
from ~ 2 X 10^ K inside of the star forming region to ~ 2 x 10^ K at a 10 kpc distance from 
the galaxy whereas the density drops from ~ 4 x 10~^ cm~'^ to less then 10~^ cm~^ value. 
Such proto-galactic winds should b e detected as sources of a diffuse X-ray emission, as it is 



the case in the local Universe (e.g . IChevalierl Il992l : [Strickland fc Stevens! l2000l : ISilich et al. 



20051 : [Strickland fc HeckmanI 120091 ): 



Att I rVAx(T,Z)dr, (4) 

Rst 



wher e n{r) is the atomic number density, Ax{Z, T) is the X-ray emissivity (see [Strickland fc Stevens 



20001 ). Rout niarks the distance at which the calculations where stopped, usually set to 10 
kpc. We set the lower integral limit to Rst assuming that the X-ray emission interior to 
it is completely absorbed by the accumulated gas. The model predicts a growth in the 
X-ray luminosity in the range from 0.3 keV - 8.0 keV as one considers larger SFRs. It is 
^ 4 X 10~^ LsF ~ 2 X lO^^erg s~^ and reaches ^ O.ILsf ^ x 10"^^ erg s~^ in the case 
of model 1 with a low SFR and model 4 with a high SFR, respectively. Note that the X-ray 
emission is concentrated towards the star forming region, where the density of the X-ray 
plasma reaches its maximum value and that in proto-galaxies with a high SFR a significant 
fraction of this emission may be absorbed by numerous dense proto-stellar clouds. 

However, in the case of SCUBA sources gravity may affect the outflow significantly. 
Indeed, the escape speed at the surface of scuba sources, Vesc — i^GMdyn/ Rsf) > may 
reach ~ 1000 km s~^, value which is approximately 10 times larger than in the case of young 
stellar clusters. In many cases it is larger than the sound speed in the thermalized plasma, 
and thus larger than the outflow velocity at the surface of the proto-galactic cloud. The larger 
impact of gravity on the flow for progressively more compact systems with the same mass is 
shown in Figure 2, panel (d), which compares the run of velocity for models 1 and 2 (solid 
and dashed lines, respectively). The maximum velocity is much smaller and the flow velocity 
drops significantly with distance to the proto-galactic cloud in the case of more compact star 
forming region (model 2, dashed line). Nevertheless, it ends up exceeding the escape velocity 
value at a larger distance from the proto-galaxy center, forming a supergalactic wind. 

Proto-galaxies which lie above the threshold line (models 4 and 5) radiate a large fraction 
of the energy input rate within the star-forming volume, what leads to a bimodal hydrody- 
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Fig. 2. — The distribution of hydrodynamical variables in the supergalactic wind outflow. 
The calculations were provided for a proto-galaxy with a dynamical mass of M^yn = 2 x 
10^^ Mq, Rsf = 2.5 kpc, SFR = 2 M© yr'^ and rjM = 0.5 (model 1). Panels a, b and c present 
the run of velocity, temperature and particle number density, respectively. The dotted and 
dashed line in panel a display the local sound speed and the value of Vesdi^), respectively. 
Panel (d) compares the velocity distribution in two proto-galactic winds emerging from 
sources of different size (model 1 and 2, sohd and dashed hues, respectively). 
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Fig. 3. — The velocity of the flow for proto-galactic sources evolving in the bimodal hy- 
drodynamic regime. The calculations were carried out for a proto-galactic cloud with a 
dynamical mass of Mdyn = 2 x 10^^ Mq, SFR = 1200 M© yr^^ and rju = 0.5. The left panel 
presents the results of the semi-analytic calculations for proto-galaxies with Rsf = 2.5 kpc 
and Rsf = 1-7 kpc - solid and dotted hues, respectively. In the right panel, the semi-analytic 
velocity distribution for model 4 (solid line) is compared with the numerical results (open 
circles) . 



namic solution ( ITenorio-Tagle et al.ll2007l : IWiinsch et al.ll2008l ). In this case radiative cooling 
depletes rapidly the thermal energy (and pressure) of the thermalized plasma in the densest 
central regions of the assembled galaxy, inhibiting the fast acceleration required to reach the 
sufficient speed to leave the star forming region. This prompts the stagnation radius, Rst, to 
move out of the starburst center as it is shown in Figure 3, left panel, where solid and dotted 
lines display the semi-analytic results for models 4 and 5, respectively. Inside the stagnation 
radius the stationary solution does not exist and thus one cannot apply the semi-analytic 
model to this region (ITenorio-Tagle et al.l 120071 : IWiinsch et al.ll2008l ). However, in Figure 3, 
right panel, we show the semi-analytic velocity distribution (solid line) in the case of model 
4 and compare this with the full numerical simulations (open circles) carried out with our 
Eulerian hydrodynamic code. The code includes also the gravitational pull from the matter 
located inside the star formi ng region and a modified cooling routine, whi ch allows for ex- 
tremely fast cooling regimes (ITenorio-Tagle et al.ll2007l : IWiinsch et al.ll2008l ). The numerical 
simulations show an excellent agreement with our semi-analytic results. Both the positions 
of the stagnation point and the velocity profiles outside of the stagnation radius obtained 
with the semi-analytic and ID simulations are in good agreement. Thus above the threshold 
line the matter injected by massive stars and ablated from proto-stellar clouds inside the 
stagnation volume remains bound and is re-processed into new generations of stars despite 
the large amount of energy supphed by stellar winds and supernovae explosions. At the same 
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time, the matter deposited by massive stars outside of this volume flows away from the star 
forming region as a supersonic wind. 

The impact of gravity becomes a crucial issue if the radius of the proto-galaxy is smaller 
or equal to the critical value, presented in Figure 1 by vertical lines for different values of 
rjid and Mdyn- This occurs when the sound speed at the surface of the proto-galactic cloud 
becomes smaller than on e half of the escap e velocity, and the nominator in the momentum 



equation (equation 12 in ISilich et al.l 120081 ) goes to zero at the surface of the star forming 
region. In this case the flow velocity cannot exceed the sound speed value, the stationary 
solution vanishes and the proto-galaxy does not form a supergalactic wind. This regime is 
illustrated in Figure 4, which presents the results of full numerical simulations for a proto- 
galactic cloud with RsF < Rcrit (model 6). Here the quasi-adiabatic wind solution for a 
proto-galaxy with Mdyn = 2 x 10^^ M©, Rsf = 2.5 kpc, SFR = 40 M© yr'^ and 7]id = 0.5 
was used as the initial condition for simulations. However, the time evolution was followed 
assuming the input parameters of model 6 (see Table 1). The initial wind solution transforms 
rapidly into a complex flow with a number of discontinuities and negative velocities inside 
the star-forming region. Note that the stagnation radius is in this case at ~ 1 kpc. However 
the matter deposited between this radius and the edge of the star-forming region is unable to 
produce a superwind and instead it cools down and ends up falling towards the center. Our 
open boundary condition does not allow for the accumulation of this gas and that leads after 
a readjustment period to a recurrent cycle in which some fraction of the deposited matter 
flrst flows away but then cools down and falls back towards of the star-forming region. This 
causes the compression and storage of the hot gas into a dense shell, which is driven inwards 
by gravity. The supersonic encounter of the outer gas with the dense shell results into the 
formation of a shock wave. This at later times {t > 30Myr, dotted lines in Figure 4) produces 
the parcel of hot gas infalling behind the cold shell. The shell drives at all times a sound 
wave into the hot inner zones, what results into noticeable enhancement of temperature and 
speeds up the infalling gas ahead of the shell as it is displayed by the dotted line in Figure 
4. The simulation ended up at ~ 30 Myrs when all matter located inside a computational 
domain is falling towards the center of the proto-galactic cloud. 

Thus compact proto-galaxies with Rsf < Rcrit trap the injected matter and are not 
able to form superwinds regardless of their energy output or SFR. A number of semi-analytic 
calculations have led us to infer that below the mechanical luminosity threshold line, Rcrit 
becomes slightly smaller for proto-galaxies with a given Mdyn and rjid and a decreasing SFR. 
However the calculations showed that the difference in the value of the critical radius usually 
does not exceed ~ 100 pc. Therefore we have adopted the value of Rcrit, calculated for a 
proto-galaxy with the threshold SFR as the gravitationally bound limit for all galaxies with 
the same Mdyn and rjid- The adopted critical radii, Rcrit, for galaxies with various Mdyn and 
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Fig. 4. — The gravitationally bound regime. Panels a, b and c present runs of velocity, 
density and temperature at t = 22.2 Myr, 29.5 Myr and 31.4 Myr (solid, dashed and dotted 
hues, respectively) in the case 6 from Table 1. 
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rjid are displayed in Figure 1 by thin vertical lines. Figure 5 shows how the critical radius 
depends on rjid and the dynamical mass of the system. 

20 1 ' ' ' ' \ ' ' ' ' \ ' ' ' ' \ ' ' ' ' I 




I \ \ \ \ \ \ \ \ I I I \ \ \ \ \ \ \ \ 

0.0 0.5 1.0 1.5 2.0 



Fig. 5. — The critical size {Rcrit) of star-forming regions as a function of the ablation pa- 
rameter rjid and the dynamical mass of the proto-galaxy. The superwind feedback mode is 
inhibited in proto-galaxies with Rsf < Rcrit- The matter returned by massive stars and 
ablated from star-forming regions remains buried inside the star forming volume and as it 
accumulates it should lead to further stellar generations. The calculations were provided for 
proto-galactic clouds with 0.3 < r]idle2 and dynamical masses of Mdyn — 2x 10^^ M© (solid) 
and Mdyn — S x 10^^ Mq (dotted) lines, respectively. 

Note that in very high redshift systems (i.e. just forming galaxies) the metallicity would 
be extremely low and hence radiative cooling will be substantially lowered. This could have 
potentially a large effect on the dynamics of the reinserted matter, favoring outflows. While 
this may be true for the first galaxies, observationally it is known that the metallicity of 
high redshift systems can be large in some cases, reaching values of several times solar. If 
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metallicity grows rapidly as the galaxy forms, then radiative cooling will be even stronger 
than that predicted in our calculations further favoring the retention of the reinserted ma- 
terial. The threshold SFR line in Figure 1 moves up a factor of 1.5 approximately in the 
case of the first galaxies with a metallicity Z = O.IZq and approximately 4.5 times down for 
older systems with a super solar abundance {Z = IOZq). Note that the critical radii, Rcrit, 
remain almost identical as one considers different metallicities. A time dependent solution 
accounting for the rapid change in the metallicity of the reinserted matter will be the subject 
of a forthcoming communication. 



4. Discussion 

We have shown here that the thermalization of the kinetic energy provided by vigorous 
star formation in young forming galaxies may lead to three different hydrodynamic regimes, 
depending on the rate of star formation, the proto-galaxy total mass and radius, and the rate 
of mass loading from proto-stellar clouds. Large galaxies with low SFRs and small ablation 
coefficient r]id form supersonic winds which carry from the star formation regions the matter 
returned by massive stars and that ablated from proto-stellar clouds. Similar galaxies located 
in the SFR - Rsf parameter space above the threshold line lose via a superwind a fraction 
of the deposited matter. The matter deposited by massive stars and that ablated from 
star forming regions in the inner zones of such galaxies becomes thermally unstable due to 
strong radiative coohng, accumulates and is to be re-processed there into secondary star 
formation. Finally, the thermal pressure in compact sources with radii Rsf < Rcrit is unable 
to withstand the gravitational pull of the galaxy. In such cases proto-galaxies retain all the 
reinserted and ablated matter within the proto-galaxy volume and do not form supergalactic 
winds. 

The value of the ablation parameter r]id remains free in the theory. However, one can 
get an idea about which values of rjid are reasonable considering sources without a secondary 
star formation, which evolve in the superwind regime. Then one can notice that in the case 
of star formation with a constant SFR, which terminates when the initial gas reservoir is 
completely exhausted, the global star formation efficiency, e*, defined as the ratio of the 
stellar mass, M*, to the initial mass, Mpc, of the proto-galactic cloud at this moment will 
be: 

In equation ([5]) the stellar mass, M^:, was calculated as the difference between the mass of 
stars formed during the evolutionary time t and that reinserted by supernovae explosions 
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and stellar winds: 



M, = {SFR - Msf) X t = (SFR - i ' >< 



and the initial mass of the system, Mpc, is: 



MpG = (1 + r^irf) X SFR X t. 



(6) 



(7) 



The star formation efficiency would be then = 90% if rii d = and approaches ~ 30% 
value required to form a gravitationally bound system (e.g. iGeyer fc BurkertI 120011 ) when 
rjid = 2. Note, that the required star formation efficiency may be smaller, and thus the 
upper h mit for rjid larger, if one co nsiders a slow expulsion of the injected gas from the 
system ( iBaumgardt fc Kroupal 120071 ). 



The predictions are thus that massive star-forming proto-galaxies with large star forma- 
tion rates similar to those detected in SCUBA sources (> 10^ Mq yr~^) evolve in a positive 
star-formation feedback conditions: either in the bimodal, or in the gravitationally bound 
regime. Only proto-galaxies evolving in the bimodal regime wi ll form supergalactic winds 
as it is in the case of submillimet er galaxies SMM J 14011-I-0252 (INesvadba et al.l 120071 ) and, 
probably, SMM J221726+0013 (IBower et al.ll2004j ). Inevitably then, matter accumulation 
would follow in the central zones or in the whole proto-galactic volume. Radiative cooling 
would then reduce the injected gas temperature, what would promote an even stronger cool- 
ing and recombination, making the accumulated gas an easy target of the UV radiation field. 
Photoionization of this gas is to set an equilibrium temperature (Thu < 10^ K), causing 
it to become Jeans unstable, leading unavoidably to its collapse and to the formation of 
new stars. Many stellar generations are expected in this scenario, until most of the mass, 
through its continuous recycling, has been converted into low mass stars with M < 7 Mq. 
The resultant stellar populations and the ISM would then show a large metallicity spread. 

Consequently, if the formation of large stellar spheroids (galaxy bulges or elliptical 
galaxies) occurs through a process of rapid mat ter accumulation and further conversion of 
this matter into stars (e.g. ISwinbank et al.l 120061 ) . which would imply a large star formation 
rate, the expectations are thus that little or none of the returned matter, through winds and 
SN explosions, is going to be ejected out of the system. Instead it is to be reprocessed into 
further episodes of stellar formation. This implies that the largest episodes of star formation 
would leave little trace of their stellar evolution into the intergalactic medium leading instead 
to a fast metal enrichment of the interstellar gas, as observed in high redshift quasars (e.g. 
Hamann &: FerlandlllQQQt lJuarez et al.ll2009l ). The hydrodynamics of star-bursting galaxies 
with a central supermassive black hole will be the subject of a futher communication. 
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